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We present a methodology for computing photocurrent production in molecular semiconduct- 
ing molecules. Our model combines a single-configuration interaction picture with the Schwinger- 
Keldysh non-equilibrium Greens function approach to compute the current response of a molecular 
semi-conducting wire following excitation. We give detailed analysis of the essential excitonic, 
charge-transfer, and dipole states for poly-(phenylenevinylene) chains of length 32 and 48 repeat 
units under an electric field bias and use this to develop a reduced dimensional tunneling model 
which accounts for chain-length and field-dependent behavior. 



I. INTRODUCTION 

Advances in nanotechnology have lead to the fabri- 
cation of the devices which length-scales smaller than 
the mean free path of an electron* 2 Work in this di- 
rection as pushed the scale of an individual "device" 
to the molecular scale through reported measurements 
of electronic transport through carbon nanotubes, self- 
assembled monolayers of conjugated polymers, and in- 
divisual molecules. — Moreover, there have been recent 
measurements of electrically induced light emission from 
individual single- walled carbon nanotubes^, as well as in- 
credible progress towards the fabrication and synthesis of 
direct-bandgap nanowires and super-lattices with novel 
optical-electronic properties. 

As this technology continues to press towards the 
molecular level, energy quantization, phase coherences, 
and electron-phonon coupling play increasingly impor- 
tant roles in the properties of the device. The cen- 
tral theme and challenge in designing molecular elec- 
tronic components is the manipulation and control of 
charge flow through single molecules and molecular as- 
semblages. Predicting and understanding electronic cur- 
rent flow through molecular systems from an atomistic 
and first principles point of view presents a formidable 
theoretical challenge in that in requires the extension 
of standard quantum chemical methods that are well 
suited for bound state problems to solve non-stationary 
and many-body scattering in the continuum. The diffi- 
culty arises in how to impose open boundary conditions 
and steady flow-boundary conditions in a computation- 
ally feasible way. 

By in large, the theoretical description of charge trans- 
port in a molecular device can be cast in the language 
of quantum scattering theory in both time dependent 
and time- independent forms. The physical picture is 
that given by Landaues^ in which current through a 



molecular wire is generated by ballistic charge carriers 
scattering through the molecule. Hence, conductivity is 
related to the transition matrix matrix, T. For the time- 
independent case, the formal treatment involves the use 
of non-equilibrium Keldysh-Green function treatments 
to compute the T-matrix4i§*2i Within the "scattering" 
approach, one derives a self-energy contribution to the 
molecular Hamiltonian due to the contact with the leads 
attached to the wire. This leads to a non-Hermitian 
Hamiltonian with complex eigenvalues. In contrast to the 
Landauer approach, Koso\s2*i2*ii establishes a molecular 
wire as a non-equilibrium steady-state system by append- 
ing a Lagrange multiplier to the molecular Hamiltonian 
to constrain the current passing through the system to 
be some desired value. As a consequence, the modified 
Hamiltonian remains Hcrmitian, continuity is rigorously 
enforced, and the bias is determined as the difference be- 
tween the Fermi energies of the system under forward 
and reverse current. 



In this work, we focus upon electronic transport 
through semiconducting systems where Coulomb inter- 
action and electron correlation effects become important. 
We address the issue of current due to photoexcitation 
of a conjugated polymer "wire" . We focus on the case 
where the bias applied to the semiconductor is insuffi- 
cient to produce current so that all of the charge-carriers 
in the system are due to the photo-excitation process. 
We adopt the Newn-Anderson22iSi model for a two band- 
semiconducting polymer wire and use this description to 
compute the current response of the wire. 
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FIG. 1: Scheme of the one-particle description and definition 
of injection energy with respect to the band-center. 



II. THEORETICAL MODEL 

A. Metal-Polymer-Metal Hamiltonian 

We consider here a rather idealized scenario in which 
an jV-site polymer chain bridging two metal electrodes. 
This is idealized in the generic sense because in the ma- 
jority of photo-current experiments, the polymer chro- 
mophorcs in the sample are either randomly distributed 
in the sample or have more or less random orientation 
with respect to the applied bias field. Moreover, with the 
exception of experiments on single molecules or molecules 
oriented on a surface, the molecules themselves are not 
chemically bound to the electrodes. Nonetheless, we 
recognize these as possible limitations of our theoreti- 
cal model. The formalism we give here is adapted from 
the typical single-particle perspective for the case of elec- 
tron/hole interactions. By assuming that the coupling 
to the continuum is a single-particle interaction, we can 
adopt a general scattering-theory perspective to describe 
current producing state in the presence of electron/hole 
excitations. 

To begin, we partition the state space of our system 
into three domains, Qlu and Qp a which span the states 
of the left and right electrodes and P which spans the 
electronic states of the bridging molecule and be pro- 
jected into separate Fock spaces representing different 
total numbers of elementary electron/hole excitations. 

P = P a + Pi + P 2 + ■ ■ ■ 
Thus, the full Hamiltonian has the structure 



H La V PL 
11= | V PL H P Vpr 

VpR H^a 



(1) 



where the diagonal terms are the hamiltonians for the 
uncoupled subsystems and Vpl and Vpp are the inter- 
actions between the polymer and the left and right elec- 
trodes. We will assume that number of excitations within 
the total system is conserved so that Vpl and Vpr can- 
not couple different polymer Fock spaces and are thus 
single particle (electron or hole) operators. 

Thus, within the space of single excitations, the poly- 



H P = 22(F mn + V mn )\m)(n\ 



(2) 



where \m) = |m e m^) is mono-excited electron/hole con- 
figuration, V mn is the two-particle matrix-element, and 
F mn is the matrix of the one-particle band-structure op- 
erator, / (which also includes any external field bias). In 
the absence of the electrodes, this is given by 



(mh\f\n,h) 



(3) 



When the coupling to the electrodes is introduced, these 
terms need to be modified to include the one-particle 
polymer-electrode coupling. For this we need to make a 
number of assumptions regarding the form of this cou- 
pling. First, we assume that the polymer-electrode cou- 
pling occurs only at the terminal sites of the polymer. 



V L p = rjkm e m h (\k)(l e mh\ + \k)(m e l h \). 



(4) 



Using the Feshbach method 2-1215 and assuming that the 
density of states of the electrodes is that of a 3D metal, 
we can derive reduced one-particle matrix elements, F 
for the terminal sites of the polymer. 



Fn — Fn — Si(-E) 
Fnn — Fnn — Sjv(-B) 



(5) 
(6) 



with the remaining terms left unmodified. 
Ak(E) — iKk{E) is the complex (retarded) self-energy 
contribution from the electrodes. (Note: T, K = (S^-)* 
is the advanced self-energy.) For these, we adopt the 
Newns- Anderson 20 - 21 model and write 



A K (E) = Z± 




\E/H < 1 
E/2"/ < 1 

E/2 1 > 1 



\E/2j\ < 1 
otherwise. 



(7) 



(8) 



where E is the "injection energy'' for an electron (or hole) 
to tunnel between the metal and the polymer measured 
relative to a common Fermi energy, Ep 1 47 is the band- 
width of the metal reservoir, and (3k is the chemisorption 
coupling between the polymer and electrode. This sys- 
tem is illustrated in Fig.^ For a metal/polymer contact, 
we take 47 = 40 eV as the band-width of the electrode 
and with a chemisorption coupling of /3k = 0.5eV. For 
the unbiased case, the injection energy for an electron or 
a hole is half the band-gap of the polymer, E = 2.5eV for 
PPV so that both the bonding and anti-bonding orbitals 
of the 7r-electron system fall into the band-width of the 
metal. Thus, I-B/27I < 1 and we can approximate the 
retarded self energy as 



(9) 



3 



which indicates that the coupling to the electrode should 
produce only a minor perturbation to the real energy of 
the system and that its only effect is to provide a sink for 
the e/h excitation to decay. We will also assume through- 
out that the applied bias is insufficient to cause direct 
current flow through the system. The Newns-Anderson 
model was also used by Mujica et aZ i 1 ?! 14 for the descrip- 
tion of scanning tunneling microscopy (STM) current in 
molecular imaging of one-dimensional systems to encom- 
pass the more general process of electron transfer be- 
tween two reservoirs. 

Our primary assumption throughout is that the self- 
energy contribution is a single-particle term. Thus, we 
can include the effect of the coupling to the metal leads 
into a configuration interaction theory by including this 
term in the generation of the Fock-matrix 

F mn = fe(m|(/+S)|n)-5 mn (m|(/ + S)|n) (10) 
where 

(m\Z R \n) = 5 ro „(£„i£f(e + /z) + <5 n jv£$(£-M) 
<m|£*|H) = fe(%Ef(- £ -/i) + ^iv^(-s-/i)) 

were ±e is the injection energy for an electron (+) or hole 
(-) onto the terminal site in the absence of an applied 
bias and /i is the shift of the bare site energy due to the 
applied field. For this we take the Fermi energy of the 
system as a common reference. Including the self-energy 
into mix produces a non-Hermitian complex symmetric 
CI Hamiltonian which can be diagonalized via unitary 
transformation yielding complex eigenvalues, e Q , 



(H P + J: R )\^ a ) =e a \ip. a ). 



(11) 



Since the polymer/electrode interaction is at the ter- 
minal ends and only involves single-electron terms, we 
have an exact relation between the net electron/hole pop- 
ulation for a given CI state on the terminal sites and 
complex energy shift component of the CI eigenvalue, 
£a = £a + 5e a + i\ a where the real energy shift and 
imaginary energy are related to the real and imaginary 
components of the self-energy, 



8e a 



A K(e La + e Ra + h+ a + h+ a ) 



A K(e La + e 



Ra 



h + Ra ) 



where, 



riT 

h t a = E \{ml\ip a 

m 
777 



(12) 
(13) 

(14) 
(15) 
(16) 
(17) 



eigenstates with more total electron or hole population on 
the terminal sites are more coupled to the metal contin- 
uum than eigenstates with less amplitude on the terminal 
sites. 

Once we have obtained the complex-eigenspectrum, we 
can construct the retarded Greens function for e/h scat- 
tering through the molecule. 



G R {E) = (E-H + ir])- 1 

= {E-Hp-Y, 11 )- 1 . 



(18) 



Because Hp + T, R is not a Hermitian operator, its eigen- 
vectors do not form a orthogonal set and we need to 
compute the eigenspectrum of the adjoint matrix 



(H P + Y, A )\^ a ) =t\ip a ). 



(19) 



Taken together, the eigenvectors {|^q}} and {\ip a )} form 
a bi-orthogonal set 

5~Jv'a(n)*V /J (n)=« a/ ,. (20) 



Consequently, the retarded Greens function is con- 
structed via 



G (n, n') = 



E 



(n\ip a )((p a \ri) 

E — e n 



(21) 



In a similar vein, the spectral response, which gives the 
generalized density of states inside the polymer taking 



N(E) 




2 3 4 5 6 7 



are the populations of the electron or hole on the termi- 
nal (i.e. right and left) ends of the polymer. In essence, 



FIG. 2: Density of states for PPV48 under 0.03eV/unit bias. 
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the contact with the metallic continuum into account is 
given by 



A = i{G R -G A ) 
= G R TG A 



(22) 



where V is the energy broadening taken as as the differ- 
ence between the retarded and advanced self-energies 



T = i(Z R - E A ) 



(23) 



Pulling these two together yields the usual expression 

(p a (n)ip a (n') 



(24) 



where e a is the "bare" CI eigenvalue of Hp in the absence 
of the metal contact. Lastly, the density of states is given 
by the trace of A, 



N(E) 



-E 

2tt ^ 



(E-e a + 6e a ) 2 + {\ a /2f 



(25) 



Thus, each CI state in becomes broadened by the contact 
with the metal implying that the states take on a finite 
life-time in addition to being displaced in energy. In the 
case of weak interaction, X a — > and the "Lorentzian" 
describing the spectral contribution from the a state 
collapses to a ^-function. The density of states for a 
32-repeat unit model of PPV under 0.03eV/unit bias is 
shown in Fig. [21 

In Fig. [3] we compare the distribution of the real and 
imaginary components of the singlet SCI states for a PPV 
polymer chain of 48 and 32 repeat units with a bias po- 
tential of 0.03eV/unit cell. This corresponds to a macro- 
scopic field of 1.0 x 10 5 V/ cm which is on-par with the 
fields placed across actual devices. The shading of each 
point in the figure reflects the net charge on the terminal 
ends of the molecule with white or red circles indicat- 
ing no net charge on the terminal site. Also indicated 
on this plot are a set of energies we will associate with 
the lowest exciton (XT), the lowest energy pure charge 
transfer state where the electron and hole do not reside 
on the same repeat unit (CT) , and a "dipole" state (D) 
in which the electron and hole are localized on opposite 
ends of the chain. 

We see that the energies of the lowest energy excitonic 
state are identical for both chains at 2.5eV as is the loca- 
tion of lowest pure charge-transfer state at about 2.9eV. 
Both of these sets of states are indicated by circles on 
Fig. [3J What is interesting in this regard is that there 
appears to be a progression of purely excitonic states in- 
tersecting a progression of charge-separated states. This 
progression of excitonic states are simply the center-of- 
mass translational eigenstates of the lowest energy bound 
electron/hole state. Even though the exciton states ac- 
quire quite large (in comparison) imaginary energy com- 
ponents indicating that the life-time of the states steadily 
decrease as energy increases, they remain relatively un- 
polarized by the applied field. As we shall demonstrate 
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FIG. 3: (color) Complex energy distribution for the SCI 
states of PPV32 (colored disks) and PPVas, (gray-disks) under 
0.03eV/unit bias. Note the imaginary energy component A is 
scaled by 10 3 . The shading or hue for each point indicates 
the degree of polarization of each state as measured by the 
net charge on the terminal end of the chain with open circles 
indicating no net charge on the terminal site (red for PPV32). 
The six circled energies correspond to the lowest energy dipole 
state (D), exciton (XT), and charge transfer (CT) state for 
either chain. 



next, the current produced by a given excitation is deter- 
mined by both the size of the imaginary energy compo- 
nents and by the net charge on the terminal ends of the 
molecule. Consequently, the excitonic states will con- 
tribute very little to the current where as the charge- 
transfer states will give the dominant contribution. 

In Fig. [H we plot the XT, D, and first two CT states 
for the 32 unit chain. Here, the exciton is more or less 
unperturbed by the applied field and is more or less iden- 
tical to the exciton for the polymer in the absence of the 
field. For the dipole state (D), the electron and hole are 
localized on the opposite ends of the polymer. For the 
32 unit chain, E{D) > E(XT); however, for the 48- unit 
chain the dipole state is considerably lower in energy than 
the XT. The energy of this state should scale almost lin- 
early with chain size since to lowest order approximation 
E(D) = qna as per a dipole in an electric field. The CT 
states shown in Fig. 0] are tunneling states from a bound 
CT state (i.e. the mAg) into the continuum. In fact 
this entire picture can be rationalized in a simple one- 
dimensional model of the unimolecular decay of a bound 
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FIG. 4: Electron/hole distributions for the XT, D, and two 
CT states for PPV32 at 0.03eV/unit bias. The energies for 
these states are circled on Fig. ED 



electron/hole state as illustrated in Fig. EI 



B. Current response to an excitation 



FIG. 5: One dimensional reduced model for unimolecular de- 
cay of a bound electron/hole pair in the presence of an elec- 
tric field. The total electron/hole potential is indicated by the 
thick line and the corresponding exciton, dipole, and charge- 
transfer states are superimposed. The shaded "wall" corre- 
sponds to the end of the polymer connected to the continuum. 
As the length of the polymer increases, the D state moves far- 
ther to the right and lower in energy whereas the XT and CT 
remain more or less unchanged. 



If we consider the net flow of charge from the molecule 
following promotion to some state <p a , we need to be 
cognizant of the fact that even though the lifetime of the 
state is given by 2X a /h, this represents only the net decay 
of the electron/hole excitation from the molecule into the 
metal since there is an equal influx and efflux of charge 
to and from the molecule. In other words, the current 
produced by a given state is not simply proportional to 
the rate of decay of the state into the continuum. Fur- 
thermore, in the limit that the valance and conduction 
bands have identical interaction with the metal (as in the 
case at hand) no net charging of the molecule can occur 
as a result of contact with the metal. Hence, total charge 
Q on the molecule must remain a constant of the motion. 
Since I — Q, one is tempted to conclude 1 = 0. However, 
the net current should be taken as the difference between 
its retarded and advanced components. 

I = (Q R - Q A ) = i[^ R - Z A , Q]/h (26) 
= e\T,Q]/h (27) 



Thus, the current produced at a given excitation energy 
E is given by 

1(E) = Tr([r,Q])/7i 



e 1 



^La-h + La ) + {h + Ra -e Ra )) ( 



E : 'a W^La '"La) 1 V "Ha v Ka/// oQ\ 
„ (E - e a + 5e a f + (A a /2) 2 1 ' 

For the case of electron/hole symmetry, we have the 
additional requirement 



e La h La — h Ra e Ra . 



(29) 



Thus, we can consider the current produced by a given 
state by writing 



(2e\ a /h) 



hole any w her' 



1+ on site 1 
1 anywhere 



(2eX a /h)(e- La -h+ a ) 



(30) 



This is similar to the notion by Halpern that current 
carrying states are localized at the opposite ends of the 
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FIG. 6: Current vs. excitation energy and spectral response 
for PPVis, at 0.03eV/unit bias. Here we show both the nat- 
ural line-widths of the current resulting from a single chain 
and the inhomogeneously broadened current response. In this 
case, the effect of the field is to shift the lowest current pro- 
ducing state below the lowest energy exciton so that current 
is produced without any excess photon energy above the state 
optically coupled to the ground state. 



conducting sample 2 — since these are the states with the 
most interaction with the terminal leads. For purely ex- 
citonic states, the net charge on the terminal sites will be 
exactly zero. Hence a purely excitonic species will pro- 
duce no net current. However, a charge-transfer species 
will produce current when the system is placed under a 
bias. Reversing the bias produces current in the oppo- 
site direction since the current producing charge-transfer 
states will be polarized in the reverse direction. 

In Fig|S]we show the current response for a given ex- 
citation energy as well as the optical response for PPVi% 
under 0.03eV/unit bias. From the way we have defined 
the direction of the current, a negative response indicates 
current in the forward direction, i.e. electrons flowing 
from the — to + ends and holes flowing + to — . A zero 
response indicates that the electron current exactly bal- 
ances the hole current, so no net charge flows through the 
system as a result of excitation. At roughly 5.5 eV (center 
of the single-excitation band) the current response passes 
through zero and we see current going against the direc- 
tion of the bias field as the energy increases. Since this 
is well above the optical response of the system which 
peaks at about 2.5 eV, this reverse current most likely 
would not be observed in any molecular system. 



C. Effect of Bias Field 

Finally we turn our attention towards what happens 
to the induced current as a function of applied field. In 
Fig. [7| we show the current response for a 32-unit chain 
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FIG. 7: Current vs. excitation energy for PPV32 in two dif- 
ferent applied fields (O.OleV/unit and 0.03 eV/unit). Here 
we have included a 100 xA a inhomogeneous line broadening 
to each state to mimic the effects of environmental disorder. 

under 0.01 eV/unit and 0.03eV/unit bias and the corre- 
sponding complex CI energies in Fig. [S] First we note 
that the on-set of current shifts towards higher energies 
upon going from 0.03eV/unit to 0.01 eV/unit. This is 
commensurate with the fact that as the field decreases, 
there is proportionally less Stark shifting of the dipole 
states. Consequently, at low field we have fewer poten- 
tial current producing states at lower energies compared 
to the higher field case. On the other hand, the XT 
states remain completely unperturbed by the change in 
field. At lower field, the lowest CT states have smaller 
imaginary energies than at high-field. Again, this is con- 
sistent with our reduced one dimensional model in Fig- El: 
lower-fields produce less tunneling between the CT states 
and the continuum. 



III. DISCUSSION 

In this paper we have laid a theoretical ground-work 
for the study of current producing states in molecular 
semiconductors focusing solely up current produced upon 
electronic excitation. In a subsequent works, we will ad- 
dress the issue of DC and AC conductivity in molecu- 
lar semiconductors, especially those in which half of the 
molecule is more electron rich than the other. Such model 
pn-junctions have been produced synthetically^^ and 
we have recently studied the electro-luminescence of such 
systems —& 

Needless to say, there are several lacuna in our treat- 
ment that will need to be systematically addressed. First, 
we assume that the electric field across the molecule is 
uniform. In realistic systems, this may not be the case 
and one really should self-consistently solve the Pois- 
son equation for the field taking into account the actual 
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FIG. 8: Real and imaginary SCI energies for PPV32 in two 
different applied fields. Gray and open circles: O.OIeV/unit. 
Red and hue-shaded disks: 0.03 eV/unit). 



tween the metal leads and the polymer than is accounted 
for in the present model. Secondly, we do not currently 
account for lattice excitations or distortions due to the 
applied field or the current through the system. The cur- 
rent response as we have defined it, simply gives the co- 
herent ballistic current produced at some energy E above 
the HF ground state. We are currently working to extend 
this to finite temperature by the inclusion of phonon cre- 
ation and annihilation channels. 

In spite of these shortcomings, the basic physical pic- 
ture offered by this model is quite appealing. First of all, 
it emphasizes the fact that different electron/hole states 
have widely differing capacities for producing current fol- 
lowing a given excitation. Possessing a large imaginary 
energy alone is insufficient to produce current. The state 
must also be sufficiently polarized so that there is ade- 
quate charge separation. This raises another crucial as- 
pect in trying to use this model to understand photocur- 
rent in bulk systems. That is, states with large dipole 
moments will be poorly coupled to the ground state via 
optical excitation. Consequently, direct excitation from 
the ground state to a strongly current producing state 
is highly unlikely. What is more likely is that following 
photo-excitation to some higher-lying eigenstate, phonon 
creation kicks in and induces transitions between the ini- 
tial state and lower energy electronic states including 
those capable of producing current. 



charge distribution of the molecule and the appropriate 
boundary conditions imposed by the conducting leads. 
Similarly, the model fails to account for the image charges 
between the polymer and metal contact. Both of these 
effects require more detailed dynamical interactions be- 
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